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We report the observation of charmless hadronic decays of charged B mesons to the final state 
K + K~tt + . Using a data sample of 347.5 fb" 1 collected at the T(AS) resonance with the BABAR 
detector, we observe 429 ± 43 signal events with a significance of 9.6 a. We measure the inclusive 
branching fraction B(B + -> K + K~ty + ) = [5.0 ± 0.5(stat) ± 0.5(syst)] x 10" 6 . Inspection of the 
Dalitz plot of signal candidates shows a broad structure peaking near 1.5 GeV/c 2 in the K + K~ 
invariant mass distribution. We find the direct CP asymmetry to be consistent with zero. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



B meson decays to final states with even numbers 
of strange quarks or antiquarks are suppressed in the 
Standard Model. Such decays may proceed by the b — > 
d loop (penguin) transition, or by other processes fol- 
lowed by ss production. Hadronic b — > d penguin tran- 
sitions have recently been observed [3, 0], while exam- 
ples of ss production have been seen in various B de- 
cays 0, 0, HI- Furthermore, Dalitz plot (DP) analyses 
of B+ -► K+K+K- @, 0| and B° -> K + K~K° Q 
have seen anomalous excesses of events at low K + K~ 
invariant masses, the origin of which has aroused consid- 
erable interest among theorists @, 11. 12. as it is 
of great importance in the understanding of low energy 
spectroscopy 1J] . Understanding the production mecha- 



nism of charmless B decays to such multibody final states 
is therefore a priority. 

The decay B + — > K + K~ir + and charmless quasi- 
two-body B decays resulting in this final state have 
not been previously observed. The current experimen- 
tal upper limits are B(B+ -> K+R-1T+) < 6.3 x 
10" 6 [H, B(B+ -> 07T+) < 2.4 x 10" 7 and 
B(B+ -► K+K*°{892)) < 1.1 x 10~ 6 0, all at 
90% confidence level. Such decays play an important 
role in analyses based on flavor SU(3) that can limit 
the allowed values of the deviation of sin(2/3 ) mea- 
sured in hadronic b — ► s penguin modes to the ref- 
erence value obtained in b — > ccs transitions such as 
B° — > J/tJjK® [18]. Various theoretical predictions give 



B(B+ -» 07T+) < O(10- 8 ) [12L l20|, l2l|, l22j, 123J and 
B(B+ -> A+A*°(892)) < O(10" 6 ) 0, M, M, Hi 
A recent phcnomcnological analysis gives a lower bound 
of B{B+ -► A+A*°(892)) > 0.7 x 10" 6 

We report herein the results of a search for the charm- 
less hadronic decay B + — > K + K~ir + [27] . The data used 
in this analysis, collected at the PEP-II asymmetric en- 
ergy e + e~ collider consist of an integrated luminos- 
ity of 347.5 fb -1 recorded at the T(45) resonance. In ad- 
dition, 36.6 fb -1 of data were collected 40 MeV below the 
resonance. These samples are referred to as on-resonance 
and off-resonance data, respectively. The on-resonance 
data sample contains (383.2 ± 4.2) x 10 6 BB pairs 0. 

The BABAR detector is described in detail elsewhere 
[30j |. Charged particles are detected and their mo- 
menta measured with a five-layer silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) inside a 1.5 T 
solenoidal magnet. Surrounding the DCH is a detec- 
tor of internally reflected Cherenkov radiation (DIRC), 



designed for charged particle identification (PID). En- 
ergy deposited by electrons and photons is measured by 
a CsI(Tl) crystal electromagnetic calorimeter. 

We select B + — > K + K~n + candidates by combining 
two charged kaon candidates of opposite sign with one 
charged pion candidate. Each track is required to have 
at least 12 hits in the DCH, to have a minimum transverse 
momentum of 100 McV/c, and to be consistent with hav- 
ing originated from the interaction region. Identification 
of charged pions and kaons is accomplished using energy 
loss (dE/dx) information from the SVT and DCH, and 
the Cherenkov angle and number of photons measured by 
the DIRC for tracks with momenta above 700 MeV/ c. We 
distinguish kaons from pions by applying criteria to the 
product of the likelihood ratios determined from these 
individual measurements. The efficiency for kaon selec- 
tion is approximately 80 % including geometrical accep- 
tance, while the probability of misidentification of pions 
as kaons is below 5 % up to a laboratory momentum of 
4 GeV/c. 

Continuum e + e~ — * qq (q = u,d,s,c) events are the 
dominant background. To discriminate this type of event 
from signal, we use a neural network [3l| that combines 
five variables: the ratio of the second order momentum- 
weighted Legendre polynomial moment to that of the ze- 
roth order; the absolute value of the cosine of the angle 
between the B direction and the beam (z) axis; the mag- 
nitude of the cosine of the angle between the B thrust 
axis and the z axis [all quantities calculated in the center- 
of-mass (cm.) frame]; the product of the B candidate's 
charge and the flavor of the recoiling B as reported by 
a multivariate tagging algorithm [32 ]: and the boost- 
corrected proper time difference between the decays of 
the two B mesons divided by its variance. 

In addition to the neural network output (NN out ), we 
distinguish signal from background events using two kine- 
matic variables: the difference AE between the cm. en- 
ergy of the B candidate and y/s/2, and the beam-energy 
substituted mass toes = \/ $1^ ~ P^i where y/s is the 
total cm. energy and ps is the momentum of the can- 
didate B meson in the cm. frame. We select signal 
candidates that satisfy NN out > 0.29, 5.272 < toes < 
5.286 GeV/c 2 and \ AE\ < 0.075 GeV. 

Another potentially large source of background arises 
from B decays containing charm mesons and charmo- 
nia. We veto B candidates with A + A~ invariant mass 
within —3 a, +5 a of the nominal D° mass, or with in- 
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variant mass of the K ir + system within ±4er of the 
mass of the J/ip or ip(2S) H. Here, a is 25 MeV/c 2 for 
D°, and 21 MeV/c 2 for J/ip and ^(25). The asymmet- 
ric D° veto is chosen to remove backgrounds resulting 
from 7r — * K misidentification. Charmonium contribu- 
tions arise mainly from the leptonic decays of J/tp and 
tp(2S), when one lepton is misidentified as a pion and the 
other as a kaon. 

The efficiency for signal events to pass the selection 
criteria is 22.1%, determined with a Monte Carlo (MC) 
simulation in which events uniformly populate the Dalitz 
plot. The only selection requirements that exhibit any 
strong dependency on the DP position are the track pre- 
selection (due to the reduced acceptance of low momen- 
tum tracks), and charm and charmonia vetoes. The av- 
erage number of B candidates found per selected event 
is 1.12. In events with multiple candidates we choose the 
one with the highest probability of a fit of the three tracks 
to a common vertex. In about 1 % of signal events the 
B candidate is misreconstructed due to one track being 
replaced with a track from the rest of the event. Such 
events are considered as a part of the signal component. 

We study possible residual backgrounds from BB 
events using MC simulations. We find that these can be 
conveniently divided into three categories, each having 
similar shapes in AE and toes- The first two (BB\ and 
BB 2 ) are dominated by specific decays, B + — > K + tt + tt~ 
and B + — ► K + K + K~ respectively. The third category 
(BB3) contains the remainder of the BB background, 
and is mainly combinatoric in nature. Based on our MC 
studies, the total number of BB pairs in our data sample, 
and the branching fractions listed by 3^, 34 1, we expect 
69, 255, and 528 events from the three BB background 
categories, respectively. 

In order to obtain the B + — > K + K~tt + signal yield, 
we perform an unbinned extended maximum likelihood 
fit to the candidate events using three input variables: 
m E s, AE and NN' out = 1 - arccos(2AW out - 1)- The 
NN^ ut variable is designed to allow simpler modeling of 
the strongly peaking structures near zero for continuum 
background and near one for signal. For each event cat- 
egory j (signal, continuum background, or one of the 
three BB background components), we define a proba- 
bility density function (PDF): 



Vi^VjimEs^VjiA^VjiNN, 



out J 



(1) 



where i denotes the event index. This form of the PDF 
is found to be valid since correlations among the input 
variables are small. The extended likelihood function is: 



£ = n cxp (~ nfe )ii 



(2) 



where rij is the yield belonging to the event category j. 

The signal Toes and AE shapes are parametrized with 
the sum of a Gaussian and a Crystal Ball function 35 1 



and the sum of two Gaussians, respectively. We fix 
the shape parameters to the values obtained from the 
B + — > K + K~ir + phase space MC sample. The contin- 
uum background toes shape is described by the function 
xyl — x 2 exp [— £(1 — x 2 )] , with x = 2mEs/V* and £ 
a free parameter (36[, while the continuum AE shape 
is modeled with a linear function. The toes PDFs for 
two of the BB background components are a Gaussian 
[BBi) and the sum of two Gaussians (BB2) while those 
for the AE PDFs are the sum of a Gaussian and a linear 
function (BB\ and BB2). The BB3 background cat- 
egory has the same functional forms as continuum in 
both toes and AE, and discrimination between these 
categories is provided only by NN^ ut . We use one- 
dimensional histograms to describe all NN^ ut distribu- 
tions. These are obtained from MC samples for signal 
and the BB background categories, and, for the contin- 
uum background, from a combination of off-resonance 
data and on-resonance data in a continuum-dominated 
sideband of toes and AE. 

The free parameters of our fit are the signal and con- 
tinuum yields, together with the £ parameter of the con- 
tinuum Toes shape and the slope of the continuum AE 
shape. All shape parameters and yields of the three BB 
background categories are fixed according to the MC ex- 
pectations. All NN^ ut shapes are fixed. 

We test the fitting procedure by applying it to ensem- 
bles of simulated experiments where events are drawn 
from the PDF shapes as described above for all five cat- 
egories of events. We repeat the exercise with qq events 
alone drawn from the PDF into which we embed signal 
and BB background events randomly extracted from the 
MC samples. We find negligible bias on the fitted signal 
yield in either case. 

Using the fit described above to the 16143 candidate 
events, we find 429 ± 43 signal events and 14850 ± 129 qq 
background events. The results of the fit are shown in 
Fig. Q] Both the Toes and AE distributions show clear 
signal peaks. The statistical significance of the signal 
yield, calculated from the change in negative log like- 
lihood with signal yield floated compared to that with 
signal yield fixed to zero is 12.6 a. 





5.274 5.276 5.278 5.28 5.282 5.284 

m ES (GeV/c 2 ) 



FIG. 1: Projections of candidate events onto a) toes and b) 
AE following a signal enhancing cut on the likelihood ratio 
calculated without the plotted variable. Points show the data, 
dark filled histograms show the qq background and light filled 
histograms show the BB background component. 
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We obtain the inclusive branching fraction of B + — > 
K + K~7r + using the result of the fit to calculate sig- 
nal probabilities for each candidate event [37|. These 
are divided by event-by-event efficiencies, that take the 
DP position dependence into account, and summed to 
obtain an efficiency-corrected signal yield. We further 
correct for the effect of the charm and charmonia ve- 
toes, and divide by the total number of BB events in 
the data sample. The result is B (B+ -> K+I<-it + ) = 
(5.0 ± 0.5 ± 0.5) x f 0~ 6 , where the first error is statisti- 
cal and the second is systematic. The systematic error 
arises due to uncertainties in the PDF shapes (2.8 %) in- 
cluding possible data-MC differences in the signal PDF 
shapes studied using a control sample of B + — > _D°7r + , 
D a — > K + it~; potential fit biases, dominated by the 
change in the result when the yields of the BB back- 
ground components are floated (6.1%); uncertainties in 
the efficiency, due to tracking (2.4%) and PID (4.2%); 
uncertainty in the correction due to vetoes, arising from 
the nonuniform DP structure of the signal, and estimated 
from MC simulations with different resonant contribu- 
tions (6.1 %); and the error in the number of BB pairs 
(1.1%). The significance of the signal including system- 
atic uncertainties is found to be 9.6 a from the change in 
negative log likelihood with and without the signal com- 
ponent, while varying those sources of uncertainty that 
affect the signal yield (PDF shapes and yields of BB 
background components). 

We also extract the direct CP asymmetry in the inclu- 
sive signal yield by separately fitting B~ and B + sam- 
ples. The asymmetry is obtained using Acp — jv-+l/v+ 
where N~ (A + ) is the fitted signal yield in the B~ (B + ) 
sample, corrected for efficiency and veto requirements. 
We find A C p = 0.00 ± 0.10 ± 0.03, where the first error 
is statistical and the second systematic, including uncer- 
tainties in the BB background estimation (0.02) and pos- 
sible detector asymmetry (0.02). Other possible sources 
of systematic error are found to be negligible. 

The efficiency-corrected Dalitz plot for signal decays, 
obtained using event-by-event signal probabilities, is 
shown in Fig. [2 We have checked that this technique 
correctly reconstructs the signal DP distribution using 
MC simulations in which the B + — > K + K~ir + events 
contain different structures. In the data, we see an ex- 
cess of events at low K~tt + invariant mass, and a large 
enhancement due to a broad structure at low K + K~ in- 
variant mass. To further clarify these structures, we show 
in Fig.[3]the respective invariant mass projections follow- 
ing requirements that remove low mass combinations on 
the other axis of the Dalitz plane. Approximately half of 
our signal events appear to originate from the structure 
at low K + K~ invariant mass. We have studied the Dalitz 
plot distributions of the backgrounds, which are found to 
be consistent with expectations, and do not contain any 
structures that may explain the peak in the K + K~ in- 
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FIG. 2: Efficiency-corrected Dalitz plot distribution of B + — > 
K + K~tt + decays, obtained with the /Plot technique [37| ]. 
Empty regions correspond to charm and charmonia vetoes 
while the area of the boxes is proportional to the number of 
events in that bin. 



variant mass distribution. Further interpretation of this 
structure and the rest of the B + — > K + K ~ir + Dalitz plot 
requires an amplitude analysis. 
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FIG. 3: Efficiency-corrected distributions of the B + —* 
K + K~ix + signal candidates: a) m K - 7r + with m K + K - > 
2.0GeV/c 2 andb) m K+K - with m K - n + > 1.5GeV/c 2 . These 
projection plots are obtained with the /Plot technique [37| ]. 



In summary, we have made the first measurement 
of the charmless hadronic B decay branching fraction 
B(B+ -> K+K-U+) = [5.0±0.5(stat)±0.5(syst)] x lO" 6 . 
The CP asymmetry is found to be consistent with zero. 
Inspection of the Dalitz plot of signal candidates shows a 
broad structure peaking near 1.5 GeV/c 2 in the K + K~ 
invariant mass distribution that is reminiscent of similar 
structures seen in other charmless multibody hadronic B 
decays @, 0, @, [HJ . This is likely to be of great inter- 
est for the understanding of low energy hadronic bound 
states fl3 ]. 
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